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Abstract 

We study the Aqcb/Mb corrections from subleading shape functions in inclusive B-meson de¬ 
cays. We propose a natural and smooth interpolation from the endpoint region to the full phase 
space, and derive expressions for the triple differential decay rate in B and the photon 

energy spectrum in B —> Xg'j- Our results are valid to order Aqqd/Mb for hadronic invari¬ 
ant masses of order AqqdMb and to order Aqqq/M^ for larger hadronic masses. They allow a 
systematic investigation of the transition between the separate regimes of the local and nonlocal 
expansions, and can be used to study decay distributions in any kinematic variables. We consider 
several examples of interest and point out that a combined measurement of hadronic energy and 
invariant mass provides an alternative to the extraction of \Vub\ which is largely independent of 
shape-function effects and in principle allows a higher accuracy than the combined measurement 
of leptonic and hadronic invariant masses. We perform the expansion directly in QCD light-cone 
operators, and give a discussion of the general basis of light-cone operators. Reparametrization 
invariance under the change of the light-cone direction reduces the number of independent shape 
functions. We show that differing previous results for the lepton energy spectrum obtained from 
different choices of light-cone coordinates are in agreement. 
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I. INTRODUCTION 


Inclusive decays of B mesons offer a rich environment to explore the flavor sector of the 
Standard Model and to search for New Physics in radiative decays. Moment analyses of 
B decay distributions have provided a precision measurement of the Cabibbo- 

Kobayashi-Maskawa (CKM) matrix element \Vcb\ at the two-percent level, along with an 
extraction of the 6-quark and c-quark masses, and higher order hadronic parameters QBi- 
Similarly, the study of inclusive decays B —>■ allows for an extraction of \Vub\ with 

current errors of 10%-15% 0, 

The conventional treatment of inclusive B decays relies on a local operator product 
expansion (OPE) in inverse powers of the large momentum Q transferred to the hadronic 
system 0, [ill 0, 0,0, 0], However, for b ^ u transitions, tight experimental cuts are 
needed to suppress the overwhelming charm background. They usually put the kinematics 
close to the boundary of phase space where the hnal hadronic system has large energy in the 
H-meson rest frame but small invariant mass. That is, Q lies close to the light-cone with 
being much smaller then Pb' Q- Consequently, the OPE in local operators breaks down. 
The large and small components of Q can be separated by employing light-cone coordinates. 
The local OPE can then be replaced by an OPE in nonlocal light-cone operators US , which 
only expands in inverse powers of the remaining large components of Q. 

The leading term in this so-called twist expansion has been known for some time now 
EH H, 113 • When going beyond the tree-level approximation the separation of short- 
and long-distance contributions becomes important and schematically leads to a decay rate 
of the form |2ll | 


dT = H X J ^ S, 


which is factorized into a hard contribution multiplying the convolution of a jet and soft 
contribution. The hard, jet, and soft contributions are associated with the different scales 
ml 3> AqcD^^fe 3> Aq(-;j 3 . The matrix elements parametrizing the soft contributions are 
usually referred to as shape functions. 

With the current experimental reach of precision, the investigation of subleading twist 
corrections in AgcD/'^^b has become important. They were considered at tree level for the 
photon energy spectrum in B ^ I HI and the lepton energy spectrum E and 


in B 


X^^vg. 


Baryonic decays have also been 

Xcii^i was 


hadronic invariant mass spectrum 

considered jl^. The hrst investigation of subleading twist corrections in B 
given recently in Ref. [13, and it was discovered that the matching of some subleading 
contributions in the earlier B —>■ Xuii^e result [13 are incorrect. 

Beyond the tree-level approximation, the factorization into hard, jet, and soft contribu¬ 


tions at subleading order in AgcD/’^fe was hrst worked out by Lee and Stewart [27| within 


the framework of soft collinear effective theory (SCET). They investigated the general struc¬ 
ture of subleading corrections, and gave results for decay rates to C>(Aqcd/'^ 6), including 
the full triple differential rate in R —> X„£i7i>. The latter was also derived by Bosch, Neubert, 


Liie lun i ripi e umereuLiai rate in jd —» . 

and Paz 0- An analysis similar to Ref. [27j | 


was carried out by Beneke et ah j2^, too. In 
all cases, the subleading twist corrections to the differential decay rates were still given at 
tree level, mainly due to the increased complexity of the contributions arising beyond that. 

In addition, it is usually difficult to assess how far away from their literal expansion region 
the twist or local expansions are still valid, and where the transition between them occurs. 
Having a single description for the entire phase space improves this situation. By comparing 
it with the predictions of the twist or local expansion, one can systematically investigate 
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where the corrections to the local or twist expansion results become large. 

In Ref. it was shown for the lepton energy spectrum in R —>• that the standard 

twist expansion can be modihed, such that it becomes valid over the entire phase space. In 
the present paper we extend this approach to the triple differential decay rate in R —>■ 
and the photon energy spectrum in R —> ^$1- Our results are exact to O^Kqqj^IMb) for 
hadronic invariant masses sh ~ C>(Aqcd^s) and to M^) away from this region. 

In particular, they contain the complete known result for the rate to in the 

local OPE region, plus some higher order corrections. 

The results in Refs. (271,12^ |29| were obtained by hrst matching QCD onto SCET, which 
acts at the intermediate scale /i^ = AgcoAfs and allows one to take into account perturbative 
corrections at this scale and sum logarithms between the hard and intermediate scales. In 
the second step SCET is matched onto heavy-quark effective theory (HQET). However, 
since we only work at tree level, there is no need to introduce an intermediate scale and 
go through this two-step matching procedure. Instead, we directly perform an expansion 
in QCD light-cone operators. The advantage of using QCD rather then HQET light-cone 
operators is, that it preserves the structure of the light-cone OPE, not mixing it with the 
separate expansion of QCD in HQET. It allows us to dehne shape functions in QCD, which 
automatically combine all higher order HQET shape functions that would normally arise 
from expanding QCD in HQET. 

In the following section we give the basic ingredients to our calculation and discuss 
the power counting. In Sec. IHII we discuss the general basis of light-cone operators and 
their parametrization in terms of shape functions. In Sec. IHI Cl we include a discussion of 
reparametrization invariance under the change of the light-cone direction, which reduces 
the number of independent shape functions. In particular, we show that the results for the 
B —> Xulvi lepton energy spectrum derived in Ref. [ 2 ^ and with a different choice of light- 
cone direction in Refs. j27|,|2^ are in agreement. Sec. If VI contains the matching calculation 
and the results for the light-cone OPE. The results for the differential decay distributions 
are presented and discussed in Sec. El and we conclude in Sec. EB 

II. BASIC INGREDIENTS AND POWER COUNTING 

A. Hadronic Tensor and Decay Rates 

We are interested in the semileptonic decay B —>• Xuii^e and the radiative decay B —>• 

The effective weak Hamiltonian for the semileptonic decay is 



( 1 ) 


with Pl = [1 — 75 ) 72 , from which one obtains the triple differential decay rate in the rest 


frame of the B meson 




( 2 ) 


where 


Gl\Vut\^ 

° 1927r3 


( 3 ) 
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Note, that we do not include the usual factor of ml in Fq. The momentum q = pe + is 
the total leptonic momentum, Ee and E,^ are the charged lepton and neutrino energies in 
the rest frame of the decaying B meson, and we explicitly kept all phase space limits. The 
leptonic tensor is Pl] and denotes the hadronic tensor. 

For the radiative decay B —>■ Xgj the effective weak Hamiltonian has the form 

=-^^VtbV*^Cf{mb )07 with O 7 = Y^sao,i 3 E°‘^{rhbPR + rh,PL)b. (4) 


Here, denotes the electromagnetic held strength, Pji l = (1 ± 75)/2- We restrict our 
discussion to the dipole operator O7 and neglect the tiny s-quark mass. The photon energy 
spectrum in the B rest frame is 


dT" 

dK 


= sroE^9{E^)e^e*f^W^ 


s 

a/3* 


(5) 


In this case 


0 — 


Gp I VtbVt*s ?\Cf {mb) {rUb) 

3271 ^ 


( 6 ) 


where we only included the ml from the effective weak Hamiltonian. Summing over the 
photon polarization in Eq. (jS)) yields 

The optical theorem allows one to express the hadronic tensor as the forward scat¬ 
tering matrix element 




= {tL)b, 


2M, 


B 


afSi 


(7) 


where / is either u or s. We will use the shorthand {0)b = {B{pb)\0\B{pb))/2Mb to 
denote the B expectation value of some operator O between physical H-meson states. The 
operator is dehned as the imaginary part of the time-ordered product of two effective 
weak currents. 


Ti(3 = - 


TT 




( 8 ) 


The momentum q is the momentum transferred away from the hadronic system, i.e.. 


q = P£ + Pu if = u) and q = p^ (/ = s)- (9) 

The currents for f = u,s following from Eqs. ([Q) and (ED are 

Ja = U'^aPLh, Ja = Hl^aPRh- ( 10 ) 

To write J® we already used the dehnitions from Eqs. (d and (USD below. We will skip the 
flavor label / in the following when unnecessary. 


B. Light-Cone Coordinates and Hadronic Variables 

As usual, we denote the H-meson velocity by n, and dehne light-cone coordinates by 
specifying two light-cone vectors n and h, satisfying 

77,^ = = 0, n ■ h = 2, and 2 v = nPn. (11) 


4 







We use round or square brackets on indices to denote complete symmetrization or antisym- 
metrization, e.g., 




( 12 ) 


The metric and Levi-Civita tensors are decomposed as (we use £0123 = 1) 

and Vn"l (13) 

The second relation dehnes = e^'^^^nxni^/2. A generic four-vector p can be written as 

= ^P-n^ + + Pa, (14) 

with p+ = n ■ p, p_ = n ■ p, and p^ = p^Pu- 

The direction of the light-cone is hxed by setting n = —q/|q|, i.e., 

ga = 0, ga = =t |q| and thus g_ < g+. (15) 

For f = s this means g+ = 2E^ and g_ = 0. For f = u we have g^ = g+g- and Ei, = 
(g+ -|- g _)/2 — Ei, and we choose Eg and ga as our independent variables. 

Usually, the hadronic tensor is decomposed into hve scalar structure functions. For our 
purposes it will be most convenient to decompose it according to its light-cone structure, 

= ~{V + ie)fWi - ^(p - ie)fW2 - + n^n^W^. (16) 

The structure functions Wi are scalar functions of ga and g_. In terms of these and changing 
variables from g^ and Ey to ga the triple differential rate 0 takes the form 


d3F« 

dF^^dgadg- 


(9+9- (i-W’; + - 2 aa (q+q-w^ + qlw; + qiwp) 

x%_)9(2B,-,_)%+-2£f), 


where we dehned ga = Q’lt ~ 2F^£. Integrating over E^, the double differential rate becomes 
= 8r„“(9+ - 9-)"(9+9-(W'i” + W'j” + W^) + qlw: + qlW^)e(q_)e(q^ - (18) 


dgadg_ 

For B Xg'j, the photon spectrum takes the form 


— = + Wi + 2Wi). (19) 

Usually, the hadronic tensor is computed in terms of partonic variables. To express the 
decay rates in terms of hadronic variables, the total parton momentum rrtbV — q is reexpressed 
in terms of the total hadronic momentum P = Mbv — q. For example, the light-cone 
component rrib — q+ is shifted to P+ = Mb — q+ = mb + X — q+, where A = Mb — mt denotes 
the difference between the physical i?-meson and 6-quark masses. Since A ~ C>(Aqcd), this 
change of variables yields an additional source of power corrections, which has to be taken 
care of when working to subleading order. 
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We follow a different approach and directly incorporate the hadronic variables in the 
OPE, because as discussed in Sec. HTdI they are better suited for an exact treatment of the 
phase space in the twist expansion. Usually, the OPE is constructed by splitting the &-quark 
momentum as pb = rrtbV + kb^ and expanding in kb- Instead, we use 

Pb = Mbv — Xv + kb = Mbv + k with k = pb — Mbv = kb — An, (20) 

and expand in k. That is, we shift the residual momentum k by Xv compared to the 
conventional choice kb = Pb — rUbV. This is allowed, because k is only dehned up to (P(Aqcd), 
and corresponds to constructing HQET with a residual mass term 5m = A [s^. On the 
operator level kb and k turn into 

iDb = iD — mbV, iD = iD — Mbv = iDb — Xv, (21) 

where iD is the full QCD covariant derivative corresponding to pb- 
The momentum Q transferred to the hadronic system becomes 

Q = Pb — g = Mbv — q + k = P + k with P = Mbv — q. (22) 

In light-cone coordinates, 

P± = Mb — q± = P |P|, Pl = 0, and k± = kb± — A, fc_i_ = kb_L- (23) 

The decay rates will now explicitly contain only Mb, while all dependence on mb or A enters 
as higher order corrections through k±, i.e., through the shape functions. In the local OPE 
this corresponds to reexpanding mb as Mb — X, which normally does not yield a very good 
approximation to the local result, since it introduces sizable 1/mb corrections, which are 
otherwise absent. However, this is not an issue in our case, since we are going to treat 
the complete dependence exactly, as described below. Concerning /c_, the contributions 
proportional to k'^ vanish at tree level for a massless hnal-state quark, the hrst nonzero term 
involving k- being of the form kM±_ = {kb- — X)k±. Therefore, when expanded in the local 
OPE, we effectively only expand the mb dependence of a certain class of higher order 1/ml 
corrections, which should yield a very good approximation. 


C. Definition of Power Counting 

For the purpose of our discussion we formulate both local and twist expansion in terms 
of hadronic variables. In light-cone coordinates, the local OPE is obtained by writing 

„ A k-\n>^ T-, k+\n>^ 2 n n A k+k- + k\ 

= F_ IH--h1 + —- Vk/, Q= P+P- 1 -F — H-h —^^ 

^ V P-/2 V P+/2 ^ \ P+ P- P+P- 

and expanding in powers of A = Agco/Afs, where P± and k are treated as 

P± ~ 0{Mb), k±, k^ ~ (P(Aqcd). (24) 


^ We use the label b to distinguish the conventional definition kb = Pb — rabV from ours. 


6 





We will refer to this as “local power counting”. The components of k are always CI(Aqcd), 
but the size of P± varies over the phase space. When P+ becomes C>(Aqcd), k+/P+ rsj 0{1) 
is not a valid expansion parameter anymore, and the local OPE breaks down. 

The twist expansion avoids this breakdown by not expanding the dependence of Q in 
k+/P+. Usually, the formal way to achieve this is to assign the power counting 

P_ ~ 0{Mb), P+ ~ O(Aqcd), fc±, k^ ~ (P(Aqcd), (25) 

which we refer to as “twist power counting”. As P+ is explicitly counted as (^(Aqcd), an 
expansion in fc+/P+ is forbidden. However, at the same time the validity of the expansion 
is restricted to the phase space region where P+ is small, which is called the shape-function 
region. In particular, the strict application of Eq. (123) leads to an expansion in powers of 
P+/P_, including leptonic tensor and phase space, which introduces sizable errors due to 
neglected higher order terms. 

However, we can choose a different approach, such that the twist expansion becomes valid 
over the entire phase space. The basic idea is to treat P+ as an exact kinematic variable, 
i.e., to not count it as (P(Aqcd)- At the same time we still do not expand in to avoid the 
breakdown of the local OPE. In other words, we only expand in k_ and k± from the very 
beginning. To formalize this approach we dehne the power counting 

P ~ 0(E), ^ ~ 0(E). (26a) 

Here, £ is meant to be a formal expansion parameter that counts powers of k- and k±. When 
expanding in e, we treat all other quantities, including and P+, as exact. In particular, we 
do not expand in P+/P_, as is done in the standard twist expansion. This modihcation of the 
usual twist expansion was applied in Ref. to the lepton energy spectrum in P —> Xu,c^i^e, 
where the energy release Mb — plays the role of P+. 

On the operator level Eq. (I2baj) turns into 

IP ^ 0(e), P^Oie), (26b) 

where £ now counts the number of explicit covariant derivatives of a given light-cone operator. 
This already implies that expanding to automatically contains the full result to 0{M) 

in the local power counting. We perform the light-cone OPE to 0{e‘^), that is, we obtain the 
full OPE coefficients of any appearing operator with up to two explicit covariant derivatives, 
which includes all corrections of subleading order in the twist power counting. 

The size of an actual term in the expansion, for instance fc+/P+P_, still depends on the 
region of phase space, i.e. the size of P±. Since we do not count P± in any way, powers of 
£ do not correspond to powers of A = Aqqd/Mb, which is why we use £ rather then A to 
dehne the power counting. It also means that the accuracy in A of our expansion varies over 
the phase space. 

The phase space regions where the standard local and twist expansions are valid are 
pictured in Fig.^ For illustration, we take A = 0.1 and show the regions where the respective 
expansion parameters are less then a/A. Region I with P+ ~ 0{Mb) is the region of the 
standard local OPE, and the light (green) hatched area corresponds to Aqcd/P+ < \/A- 
Region H is the domain of the standard twist expansion P+ (P(Aqcd), P+/P- ~ 0(A), 
and the light (orange) hlled area shows P+/P_ < a/A. The dark (violet) edge is the resonance 
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FIG. 1: (color online) Phase space regions in the P± plane as discussed in the text. The solid line 
shows sh/M'^ = A and the dashed one sh = The light (orange) filled region is P+/P_ < \/A, 
the light (green) and dark (violet) hatched regions are P+/Mb > \/A and P^/Mb < \/A, and the 
dark (violet) filled region is the resonance region with P-/Mb < A. We take A = 0.1 in all cases. 


region P_/Mb < A, where the inclusive treatment is invalid, and the expansions necessarily 
break down. The dark (violet) hatched area Aqcd/.P- > \/A shows the transition into the 
resonance region. In region III, the vicinity of Pj^/Mb P./Mb ~ VA, the local OPE 
is also applicable, except that the expansion is only in powers of a/A. The expansion in 
e is valid anywhere away from the resonance region, and therefore provides a natural and 
smooth interpolation between the separate regimes of standard local and twist expansion. 

To investigate the accuracy of our expansion, we write as 

Q'‘ = ip_(n'‘ + A^H‘‘ + A„»- + 2A^. (27) 


We can see that itself contains an (P(A°) piece proportional to n. Taking the square. 


Q" 


p2 



P+ + 
P_ 


P+ + k+ 

p1 'K 



{P++k+)p. 


1 + 


k^ 


P_ (P+ + M^- 


, (28) 


the leading term = 0 vanishes, and the next largest term is (P+ + fc+)/P_. The scaling in 
A for the various terms in Eqs. (127j) and in regions I, II, and III is summarized in TableU 
It shows that an expansion to 0{e‘^) is exact to 0{ls?) in region I, i.e., for sh ~ 
and to 0{K) in regions II and III, i.e., for sh ~ (P(AQCDAfs). In particular, it includes all 
standard twist corrections of (P(A), as well as the complete local (P(A^) result. The largest 
corrections occur in region III, where they are only suppressed by powers of a/A. In regions 
I and II the higher order corrections are suppressed by A. 

To explicitly see the difference to the standard twist expansion, we take a closer look at 
Eqs. (|77|l and (j2Hl)- The and k^ terms in are both twist C>(A), but 0{e) and 
In itself, e.g. when multiplied by 7 ^, the k_ and kj_ terms are either twist 0{K) or 0{e). 
These are the only terms, which have a power of £. Therefore, expanding to 0(£^) includes 
all corrections of subleading twist. In addition it includes the 0{h?‘) twist terms fc^/P^ and 
k-k±/P‘^. These are precisely the 0{A‘^) twist contributions whose local expansions contain 
a local 0{h?) term. Their inclusion achieves the accuracy to local (9(A^). Note that we do 
not claim to include all (9(A^) twist contributions, which would require to include the terms 
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region 

standard OPE sh/M"^ Pj^/Mb P-/Mb 

k+ 

p+ 

P++k+ 

k-,k± 

fc2 

P 

P 

(P++fc+)P 

I 

local 

1 1 1 

A 

1 

A 

A2 

II 

twist 

A A 1 

1 

A 

A 

A 

III 

Vlocal 

A \/A \/A 

^/A 

1 

^/A 

A 

assigned power in e: 

1 

1 

e 

£2 


TABLE I: Scaling of expansion parameters in the construction of the OPE for different regions of 
phase space. The regions are shown in Eig. ^and discussed in the text. 

k-k^, and k±^k\. But the expansion to O^e"^) is correct to 0{ls?) in region I, and 
hence to 0{K) in region III. 

The second type of terms are those proportional to P+/P^, which are twist 0{A), but local 
(9(1). Expanding these “kinematic” twist terms restricts the standard twist expansion to 
small P+. In Q^, the standard twist expansion also expands the term (P+ + /c+)/P_ (9(A). 

In contrast, we do not assign a power counting to it and treat it exactly. As mentioned 
before, all kinematic factors from leptonic tensor and phase space are usually also treated as 
kinematic twist terms and expanded in P+/P_. Since they are unrelated to the OPE, we can 
treat them exactly, too. In summary, our expansion keeps all kinematic twist contributions, 
which makes it valid over the entire phase space. 


D. Phase Space and Hadronic vs Partonic Variables 

In this section, we point out a subtlety in the treatment of the phase space. Fundamen¬ 
tally, the hadronic tensor itself only contains the overall momentum conservation, while the 
remaining phase space limits are contained as 6 functions in the decay rates. Taking the 
double differential rate m as example, these are 9{mb —p-)9{p- —p+), where p = rribV — q 
is the total parton momentum. 

As dehned in Eqs. dZD and (jHl), the hadronic tensor W-^ has support for positive and 
negative values of P_, where the negative values correspond to different physical processes. 
In particular, its support is a priori not restricted to 0 < p_. Rather, when evaluating the 
imaginary part in one has to pick out the cut corresponding to 0 < P_. In the local 
OPE the hadronic tensor contains at tree level the partonic momentum conservation 

HP) = Hp+p-) = jAppH = (29) 

\p-\ P- 

The two terms in the last expression correspond to the two different cuts. The 6 function 
sets = 0 , with which the 9 functions become 9{mb — p-)9{p_). Therefore, the hadronic 
tensor is only evaluated for 0 < p_ < mi,, which automatically picks out the correct cut. 
Using the partonic variable p in the twist expansion, the momentum conservation will be 

5{{p+-uj)p-) = ^^^~'^ ^~\ {p+-uj), ( 30 ) 

where uj is the argument of the shape functions with support —A < uj. With = a; the 
phase space limits are a; < p_ < mb, and p_ can become negative. To pick the hrst term in 
Eq. (p^ we must require 0 < p_ by hand. This results in different limits on p_ depending on 
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the sign of a;, namely, 0 < a; < p_ < mi, or a; < 0 < < mi,, which is rather cumbersome. 

One could argue that this is irrelevant, since the twist expansion is only valid for large 
p_ ~ 0{mb) anyway. Also, upon integration over p_, the difference between uj < p- and 
0 < P- will be of higher order. However, since 0 < p_ restricts P_ to A < P_, this seems to 
disallow an exact treatment of the phase space, which is what we aim for. 

The problem is that we want hadronic phase space boundaries, while the partonic vari¬ 
ables force us to expand around partonic phase space. In hadronic variables, the 9 functions 
are 9{Mb — P-)9{P- — P+). Using P instead of p, the momentum conservation will be 

y(p+- (31) 

giving the limits a; < P_ < Mb- Since the support of the shape functions is now 0 < a;, this 
again picks out the correct cut by itself. In particular, P_ can now extend into the region 
0 < P_ < A. To summarize, we obtain 0 < P+ < P_ < Mb- These are the physical phase 
space boundaries, provided we neglect the mass of the tt meson, which is much smaller then 
^QCD- At present there is no way to consistently include the effects of m.„, because the twist 
expansion can only account for the nonperturbative effects due to the initial B meson. 

III. QCD LIGHT-CONE OPERATORS 
A. General Operator Basis 

All light-cone operators to 0{e‘^) can be derived from the three kernels 
/Cj)(cu) = b6{iV++ uj)Tb, 

/C[^(cui, CU 2 ) = M(iP+ + o;i)iP^5(iP+ + uj2)Tb, (32) 

^ 2 ^ (n^i, 1 ^ 2 ; 1 ^ 3 ) = b6{iP^ -f-co’i)iP^5(iP_i_ -t-co’ 2 )iP^^(iP+ +i^ 3 )r 6 , 

where T is some generic Dirac structure, and the 6-quark fields are full QCD fields. When 
DarametrizinEr the ooerators in Sec.lIIIBIwe take iP = iD — Muv. accordinsr to Eo. (EH), but 
the general discussion in this section is independent of the specific dehnition of P and uj- 
Since fC\{uii,ui 2 ) and /C 2 (i^i,a; 2 ,a; 3 ) depend on more then one variable, their parametriza- 
tions yield shape functions of two and three variables. In Refs. 1^ such parametrizations 
are given for the SCET equivalents of these kernels. At tree level only bi-local operators 
appear, and therefores only integrals of the above operator kernels are needed. In the same 
way, the shape functions appearing at tree level only depend on one variable, and are effective 
combinations of the multivariable functions |27L l29j|. 

To save some writing, we abbreviate the Wilson lines as 6+(a;) = 6(iP+ + a;). The 
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complete set of bi-local operators to is 

Olioj) =b6+{uj)rb, 

= 6[^{u;)b6+{u;i)iV>^6+{u;2)rb, 

= S'; 23 {uj)b 6 +{ui)iV>^S+{u 2 )iV’^b+{u; 3 )rb, 

OlHiu;) = h{iV^6^{u;) ± 5+(o;)iP'^)r6, 

= 6[2{uj)b6+{uJi)iV^^iV''6+{u2)rb, 

OZ"{uj) = ^6[^{u^iVf^6+{ui)iV^6+{u2) ± 6+{ui)iV^6+{u2)iVf^)Tb, 
= biV>^6+{uj)iV'^rb, 

Ol%{oj) = H{lV^^lD'^5+{oJ)±5+{oJ)lD^lD^^T)b, 


where the upper and lower sign belongs to the hrst and second label, respectively. Note 
the particnlar assignment of the Lorentz indices for and Clg^Q(a;), which tnrns ont 

to be usefnl for parametrizing them. For later convenience we dehne Clf(a;) = bOi{ui)Vb. 
That is, we drop the label F when referring to the derivative structure only, e.g., Offo;) = 
5[^8+{uji)iV>^8+{uj2). 

The 5-fnnction factors 5^2 ('^) ^123 ('^) dehned as 


<^12 (^) — 
^123 (^) = 


dcuida;2 
■ du;35^2('^ 


8{u — uji) — 8{uj 


UJi — UJ 2 
duji8'^^{uj) 


OOl - UJ 3 


J V ^^02 1^01 / 

= [ d^,du,,d^3 ( ^ ^ , 

J V<^02<^03 <^011^03 ^01^02/ 


(34) 


with ujij = ooi — (jj and ooo = 00 . They are completely symmetric in the a;* and inclnde 
implicit integrations over a;i,a ;2 and a;i,a; 2 ,a; 3 , respectively. They satisfy 

8[,{u;)8{u;i)8{u; 2) = -8'{u;), 8'l,,{u;)8{u;i)8{u2)8{u;3) = \8"{uj). (35) 

The factors in brackets on the right-hand sides of Eqs. (El arise as the imaginary parts 


1 1 5(<d2oi) , <^(<^ 02 ) 

-fm 7 -rwT-TT =- 1 -> 

TT (<^01 + ) (<^02 + ) <^02 <^01 

_1__ <^(<^oi) <^(<^ 02 ) <^(<^ 03 ) 

TT (cUol -|- F)(a;o2 + F)(<d203 + i^) <^02<^03 <^011^03 <^01<^02 


(36a) 

(36b) 


In Refs, jl^, the right-hand side of Eq. contains an additional piece 

—7r^5(a;oi)5(a;o2)^(|^03)) which we think should not be there. Eqs. (IHF)|l are dehned upon 
integration over uj. Taking the imaginary part together with the E prescription picks ont 
the poles at cu = a;i,a; 2 ,a; 3 , and the replacements in Eqs. El are a formal way of achieving 
the same. Taking the limit ui 2 )^z on both sides, Eqs. El rednce to the n = 1 , 2 cases 

of the standard formula 
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The operators in Eq. o are not completely independent with respect to their Lorentz 
structure. An operator of reduces to one of when any derivative next to a 

Wilson line is contracted with For example, using Eq. (ESI)) we have 

n^O\^{uj) = {uOl{uj))\ n^Ol^{uj) = -ujOl{uj), = (37) 

where the prime denotes the derivative with respect to uj. This simply means that only 
and "Dx (or equivalently v ■ V and T>_l) in the operators are independent structures, which 
reduces the number of shape functions needed to parametrize the operators. The full set of 
such relations is given in the appendix. 

Another comment concerns the twist order of the operators in Eq. El- Formally, the 
factors of and ^i 23 ('^) reduce the twist order of an operators. For instance, 

are formally of leading twist. Nevertheless, the discussion in Sec. Ill (Jl shows that they do 
describe sub- or subsubleading twist corrections, because they contain explicit derivatives. 
Therefore, they must have coefficients of higher twist order. We will see an example of this 
in Sec. IIV A II This involves the standard twist power counting, for one has to consider the 
operators and their coefficients. In this respect, our power counting is more transparent. 


B. Shape Functions 

In this subsection we take iD = iD — Mbv. The parametrization of the operators also 
depends on the specific Dirac structure T. We need T = 7 " and T = 7 " 75 , and define 

0“(^) = bO,{u;)rb, Vt{uj) = bO,{u;)ri5b. (38) 

We follow the notation of Ref. [ 2 ^ where possible.^ Schematically, 

{Onco))B = (Fi, G,)iu;)v^ + (K,, M,){u;){n - n)“ L,{u;)vl, . . 

{Vn^))B = H,{u)in - vr + N,{uj)v^ + ^ > 

where (Fi,Gi)(a;) stands for Fj(a;) or Giioj). In the heavy-quark limit, the 7 " in O^iuj) 
is parallel to n". Therefore, Fiiuj) and Gi{ui) contain the leading contribution, while the 
functions Ki{uj), Lj(a;), and Mi{uj) are suppressed by I/nib because v ■ {n — v) = v ■ ri± = 0. 
They contain all higher order corrections in 1 /mb perpendicular to v°‘ that would arise from 
expanding the 6 -quark field. Similarly, the axial vector 7“75 in Vf{uj) is perpendicular to 
at leading order in l/m;,, and all contributions parallel to are suppressed by l/mb- 
Similarly, there will be 1/m;, suppressions from the HQET equations of motion, see below. 


1. The Leading Operator 

The B expectation value of the leading operator is 

mu;))B = Fo(o.)n“ -F Koiu)in - n)“, {F^{u))b = 0, (40) 

which defines the QCD shape functions Fq{uj) and Kq{uj). This is exact, i.e., there are 
no higher order corrections on the right-hand side. The support of the shape functions is 


In Ref. 


0 


the operators are defined in terms of iDb, which shifts the shape functions’ argument by A. 
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0 < a; < Mb- Strictly speaking, the upper limit is Mb rather then oo, because Eq. (PUj) 
contains no reference to the heavy-quark limit. The matrix element of Vq{uj) vanishes by 
parity invariance. 

Using iD = iVh — Xv instead of iVf, in the leading operator only shifts the argument of 
the shape functions by A, such that Fo(a;-|-A), Kq{uj + X) correspond to the functions defined 
in Ref. 2^|. Their expansion into the usual HQET shape functions 1^ is 


^o(w) = /(w-A) + ^t(a;-A) + - • 


Mu) : 


(u 


mb 


^f{u-X) +—hi{u-X) + - 

mb 


• , (41) 


which explicitly shows the l/m^ suppression of Kq{uj). The QCD shape functions auto¬ 
matically contain the appropriate combinations of HQET shape functions that arise from 
expanding the QCD fields and states. This is in fact very similar to the local expansion, 
where the parameters /i^ q are defined using the full QCD states, and thus differ from the 
HQET parameters Ai ,2 by 1/mb corrections. 

To constrain the form of Fq{u) and Kq{u) we can also directly parametrize their moments 
in HQET. Using the abbreviations 


H — ^ + 372, T2 — Ts/S -|- 74, Ao — Ai -|- Ti -|- 3(A2 + T2), Po — Pi + 3p2, (42) 


where Ai, 2 , Pi, 2 , and 7i_4 are the usual HQET parameters, we hnd 

F„(w) = - A) - - A) - - A) - - A) + ■ • • , 

2mb 6 18 

IV 1 1 ‘^Xo-po/mb ,, . po p,, .. 

M^) = - 7 - ^ {uj-X) + -—6 (w - A) + ■ ■ ■ . 

bmb bmb 

It is convenient to expand with respect to a; — A, otherwise A explicitly appears in the 
moments. This is where the 6-quark mass reappears. The normalizations of Fq{u!) and 
Kq{uj) are hxed by 6-quark number conservation, while all other moments in Eq. il receive 
higher order corrections starting at order divided by an appropriate power of mb- 


2- Subleading Operators 

For the 0{e) operators, parity invariance implies that 

{0/^{n,v))B = {OMnp^'^p))B, =-{VMnp^'^p))B, (44a) 

where we temporarily suppressed the u dependence, but explicitly showed the dependence 
on the vectors n and v- The transformed vectors satisfy rip = Vp = v^- Hence, the 
must be proportional to because s'/// = —while the must not 

contain e°///‘- Time-reversal invariance requires 

(Pi%PiP))b = (Ol,3a^(np,np))5 = ((7i,3a/,(?7.p, r>p))p 

(44b) 

{02^{n,v))B = {04afj.{np,Vp))B = —{04afj.inp,Vp))B, 

and identical relations hold for the VM{ui). The complex conjugation reverses the order of 
all derivatives in the operators, yielding the additional minus sign for and V2^{u)- 

Eqs. fj44|l show that the matrix elements of Vi^{uj) and 02'^{uj) have to vanish. 
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The remaining nonvanishing matrix elements are parametrized as^ 


(C>“^(a;))s = -[(a;Fo(a;)n" + ujKo{uj){n - n)"]'(n - v)^ 

- [(Fi - AFo)(o;)n“ + (K, - AFo)(o;)(n - n)“] V - 

{0'^^{u))b = [{uFo{u)v°‘ + uKoiu){n - v)'^]{n - vY ^ 45 ) 

+ [(F3 - AFo)(a;)n“ + (F3 - AFo)(a;)(n - , 

and we already took into account the constraints from Eq. The (5 'i2 (<^) inside C>“^(a;) 

makes it formally twist 0{1). Mainly for cosmetical reasons, we want its shape functions 
to be of the twist order at which they actually appear, which is why we use derivatives of 
shape functions to parametrize the operator. Eqs. (03) are chosen such that 


{v,07{u;))b = -(Fi - AFo)'(o;)u“ - (Fi - AFo)'(o;)(n - u)“, 

{v^OT{uj))b = (F 3 - AFo)(o;)u“ + (F 3 - AFo)(o;)(n - u)“, (46) 

= Ls^uj), 

{i£±a,V7iu;))B = 

To leading order in l/rufe, ^ 4 ( 0 ;) equals —hi{uj — A) of Ref. Hi- The Fi^ 3 (a;) and Li, 3 (ci;) 
are suppressed by l/rufe, as argued before. In the heavy-quark limit the HQET equations of 
motion imply F^{ijj) = 0, and therefore F^{ui) is also suppressed by l/rub, which is why we 
choose {n — v)^ and as independent vectors in Eqs. 03). 

Considering the 0{e^) operators, by the same arguments as in Eqs. (14411 . the only non¬ 
vanishing matrix elements are 


{oi‘:^(‘^))B, (KW) 


B- 


(47) 


Due to the three indices the decompositions become rather lengthy, so we will not write 
them out explicitly.^ Instead, we define the shape functions as in Eqs. (jlF)|l by projecting 
out the independent Lorentz structures, which is done in the appendix. Here, we only list 
those needed in the following, 

{vaV±,^0^‘'^''\u;))B = - ALi)"(a;), 

{VaV±t^uO'^^^’'\uj))B = 6 * 5 ( 0 ;), {{n - v)aie±f,^V^'^^''\u}))B = ^ 5 ( 0 ;), 

(h^D±a.PiT(a;))B = (Fio - AF4)(a;). 

1-1 

The functions ( 65 ( 0 ;), H^{uj) equal 6 * 3 ( 0 ; — A), H^{uj — A) of Ref. [ 2 g^, and to leading order 
in l/rub, 6 * 2 ( 0 ; — A), ^ 2 ( 0 ; — A) of Ref. jl^, respectively. 


) there corresponds to 2 i? 4 (a;) here. 

® The numbering of the operators is changed to account for 03 ^ 4 ( 0 ;). The operators 03 ( 0 ;), 7 ^ 4 ( 0 ;) of Ref. 
correspond to 05 ( 0 ), V^iuj) here. 


^ The definition of ^ 4 ( 0 ;) is slightly different in Ref. 0 , R. 
^ A complete decomposition for is given in Ref. 


0 
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C. Reparametrization Invariance 


The operators C>o, 3 , 4 , 5 (<^) correspond to the operator basis originally introduced in 
Ref. and appear in the OPE for the triple differential rate, see Eq. dnn) below. The 
operators (Po,i, 2 , 5 (<^) are the complete set of o per ators needed for the direct computation 
of the lepton energy spectrum in i? —> Xuii'e j2g. There, the momentum of the charged 
lepton is used to dehne the light-cone direction, i.e., one chooses pi = Ein. For B Xgj 
with q = p-y this choice is equivalent to ours in Eq. For the triple differential rate with 
q = Pe + Pu the two choices of light-cone directions are rotated with respect to each other 
by an angle depending on the three-momenta of charged lepton and neutrino. 

The lepton energy spectrum is independent of the choice of the light-cone direction. 
Therefore, integrating the triple differential rate should give the same result as the direct 
computation. Since the two approaches require different subsets of operators, not all op¬ 
erators in Eq. (El can contain independent nonperturbative information, and therefore, 
some shape functions appearing in their parametrizations should be related, beyond simple 
relations like Eq. El- The shape functions basically describe the momentum distribution 
of the B meson, and because it has zero spin, they cannot contain any spatial directional 
information. However, the light-cone coordinates separate the spatial directions into ri± and 
n — n. Hence, shape functions referring to r/j, and n — n should be somehow related. 

The independence of physical quantities on the choice of the light-cone direction is de¬ 
scribed by reparametrization invariance (RPI). More generally, there are two types of am¬ 
biguities related to RPI in our setting. First, the ambiguity in the decomposition of the 
heavy-quark momentum leads to the well-known reparametrization invariance of HQET 
Its implications for the twist expansion have been studied in HQET in Ref. [s^. The 
authors there consider the case of a single outgoing particle with g_ = = 0. Their results 

are thus not applicable to the triple differential rate, and it is not surprising that they do 
not hold in this case. 

The second type of ambiguity arises from the arbitrariness in the definition of the light- 
cone vectors. This has been studied in some detail in SCET 34 ^ where it places many 

constraints on the form of allowed operators. Following Ref. u , there are three classes of 
transformations that preserve the fundamental properties = 0 and n • h = 2, 


(I) 


n (5_l, 
n. 


( 11 ) 


n, 


n + 5\ 


(HI) 


(1 -F a)n, 

(1 — a)h. 


which are generated by the five infinitesimal parameters {5±, a}. The RPI transformation 

studied in Ref. corresponds to a combined action of (I) and (HI). 

We want to study the effect of rotating the light-cone direction, while keeping 2v = n + n 
fixed. Thus, we set a = 0 and and consider the infinitesimal rotation 5r 

n ^ n + Sru, n ^ n + Sru with Sru = —Sru = (49a) 

under which 

5rv = 0, 5R{n -n) = 25±, 5rp^^ = {n - ^rS^^ = {n - n)^^e'l^5^x. (49b) 

The last two transformations can be found by requiring that metric and Levi-Civita tensor 
stay invariant. Similarly, any four-vector is invariant under 5r, only its light-cone compo¬ 
nents change according to Eqs. (Hi. 
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We also need the transformation of the Wilson line 6+ = + a;), 


<5R(5+(a;) — —(50) 

which can be found using 5+(a;) = (—l/ 7 r)Im(iP_|_ + a; + i£)“h Note, that need not 
be formally 0{A), but can be 0(1), because 6r only rotates the light-cone components of 
k into each other, but leaves v invariant. In particular, the transformation (EDD does not 
change the twist order of an operator. However, it connects different orders in e, which 
yields constraints on the shape functions arising from operators of different order in e. 

The Dirac structure T and the 6 -quark helds are unaffected by 6r, so Eq. (IKn|l yields 

6rOI{u;) = (51a) 

and for the 0{e) operators 

5rO\^{u;) = 5rOI%{uj) = . (51b) 

Similarly, the operators of 0(£^) are transformed into 0{e^) operators. Taking the B expec¬ 
tation values of Eqs. we can pull 5r out of the matrix elements, because it has no effect 
on the H-meson states, as well. Thus, the same relations also hold for the B expectation 
values, which reduces the number of independent shape functions. 

Taking T = 7 " and employing Eqs. iOl) and (03), the matrix element of Eq. I51al yields 

5r{Fo{u;)v^ + Ko{uj){n - n)“) = K,{uj)5l = (52) 

from which it follows that 

L((a;) = 2Kq{ijj). (53a) 

This relation has the expected form, since Kq{uj) and L[{uj) are proportional to {n — v) = 
{n — n)/2 and ri±, respectively. Similarly, the relation between (9“^(a;) and yields 

L' 2 {ijj) = 2Ki{uj) and G 2 {uj) = —2(a; — A)Fo(a;) — 2Fi(a;). (53b) 


27, 2>i[2 


The remaining relations following from Eqs. (ED are given in the appendix. 

Writing the second relation as G2{uj) + 2Fi(ad = — 2(a; — A)Fo(a;) one can easily see that 
the B —> Xuii^e lepton energy spectrum in Ref. [2g expanded to subleading twist agrees with 

The appearance of the different operator structure 

but 

arises from choosing the light-cone direction to be parallel to the lepton momentum. We 
disagree with the statement in Ref. 
results. Part of the reason why Ref, 


the results obtained in Refs, 
in Ref. is not related to the use of QCD vs HQET helds as presumed in Ref. 


that this choice can lead (by itself) to incorrect 
obtained an incorrect result is that it tried to 



match on an (for this choice) incomplete operator basis. 


IV. THE LIGHT-CONE OPE 

In this section, which is mainly technical, we compute the light-cone OPE of the current 
correlator in Eq. (jS)) to 0{e^) in the power counting (1^ . We start with considering 
generic currents J“ = 6P"/ and = fT^b with two arbitrary Dirac structures P" and 
P^. For completeness. Sec. IIV Al contains the actual matching calculation. In Secs. IIVBI 
and IIV (1 we give the results for general currents, Eq. and semileptonic and radiative 
currents, Eqs. dZOD. The latter are used in Sec JVAl to obtain the hadronic tensor. 
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FIG. 2: Tree-level Feynman diagrams for the zero- and one-gluon matrix element. 


A. Matching Calculation 

1. Zero-Gluon Matrix Element 


The matrix element of between &-quark states with momentum pb = Mbv -|- A; is 
shown on the left in Fig. |21 With P = Mgv — g it is 

(6|T“^|6) = -ilm (^UbT^^^-L^Yf^Ub^ = UbAo^.T^^r^/^Ub- (54) 

Using the shorthand + k^), the expansion of to O^e"^) is 

= n^5u + ^ (^n% - + 0{e^). (55) 

The derivative is with respect to the argument of the 6 function. There are no terms 
proportional to fc_ and /c^. By hrst expanding the 6 function in only, one can see that 
there are no contributions proportional to kA. 

Since Oo{uj) is the only operator at order we can extract its coefficient from Eq. (ESI) 

= j du;6{P+-u;) — {b\Oo{u;)\b) + 0{e), (56) 

where {b\OQ{uj)\b) = (5(A;+ -fa;). At higher orders, this extraction becomes ambiguous. For 
example, by partial integration we can rewrite kA^b'^ as 


kl5\P+ + k+) = j duj6\P+ - u;)kl6{k+ + ^) = j dLj6{P+ - uj)kl6\k+ + u) 
= [ duj9{P+ - uj)kl6"{k+ + uj). 


(57) 


The first expression requires 08 , 9 ( 0 ;), the second 05 , 0 ( 0 ;), and the last 02 ( 0 ;). The operators 
themselves correspond to twist 0{h?) through 0(1). The difference in their order is canceled 
by their coefficients in the convolution, which are of relative 0(A“^) through 0(A), so that 
the total order is the same. In our case the convolution always involves 5(P+ — o;). 

Because k and S{k+-\-uj) commute, while iV and 5{iD^Puj) do not, the zero-gluon matrix 
element cannot distinguish these operators, and thus only fixes a linear combination of their 
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Wilson coefficients. This is avoided by either directly expanding the propagator l/(.^ + i^), 
or gets resolved by the one-glnon matrix element, which can distinguish the operators. For 
terms containing this means that one is a priori not allowed to use the HQET equations 
of motion to replace = —fc+, because in some operators the derivatives are separated 
from the 6 -quark helds by Wilson lines. 


2. One-Gluon Matrix Element 


The one-gluon matrix element is depicted on the right in Fig. |21 It has an additional soft 
background gluon, which we take to be in the initial state with momentum I — k. Working 
in light-cone gauge = 0 , we have 


mnr,) = -hm (58) 


where 


<7. = -(f+ ^) 


^ is) 

6 [{P + k^] - 6 [{P + If] 


(59) 


{p + iy - {p + ky 

We absorbed all factors from the gluon vertex into the polarization vector e = such 

that e corresponds to a covariant derivative iD = (i9 — Mbv) -|- gA. 

According to Eq. 


we expand this in /c_, e_, and /c±, ej_. Employing the identity 

to reduce the product of three 7 matrices, we hnd to 0 {e'^) 


2^5* = -ie)'yey6k±Si) - n>^{g - ig)l'^(/ ± k)yex ^^_ ^ 

- ^ -F \eYl{eukx5k + YexSi) (61) 

-f (7 — ie)Y ((e_fci/ -|- l-eY)5k ± -l- 0(eY. 


Here, Sk = 6 (T+ -|- /c+), 61 = S{P+ -|- /+), and the upper and lower signs belong to g± and 
ie:^, respectively. 


3. Four- Quark Contributions 


We include the four-quark matrix element on the right in Fi g El which was hrst considered 
in Ref. ji^. Its size has been subject to some recent discussion [27 . Hllllal. Although the 
corresponding operator is formally of subleading twist and third local order, it is unclear at 
present how well, or if at all, this represents its actual size. We thus refrain from assigning 
it a power in e. Instead, we treat it as a separate contribution and keep only its leading 
twist term. The graph on the left in Fig. El has no imaginary part contributing to the triple 
differential rate. Its contribution to siiigle differential spectra can be computed, where the 
lepton or neutrino lines are connected il¬ 
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FIG. 3: Four-quark diagrams. On the right we only show the routing of residual momenta. 


We route the momenta such that the quarks carry residual momenta ki and k 2 , and the 
gluon residual momentum 1. To leading twist the matrix element of then becomes 


{bf\T^f^\fb) = —lm(-g^u,T 


k 2 IS 


Uf 


LLV 

VI 


{P + If + i£ 


—Uft 7 ^ 


^ P P is 


r^Uf 


(62) 


where i,j are color indices. The second form is obtained by noting that for the currents in 
Eq. (fTn|l M 4 = and hf = are left-handed spinors and employing the 

Fierz identity 

{u^l±\U2){u^7±Ui) = n(AU^)(hf7^U4)(hf7''M2)- 
The latter follows from contracting the general Fierz identity for left-handed vector currents 


{uhP)(upu^) = -(rr^- - >7>r - ifp - iE'“'"")(77A«i)(fih.« 


The expansion of to leading twist i 


IS 




Hg 


P_ 


n 


Jki 


P 




+ 


Si 


(P P k2)+{P + 0-1- (-P + 0-i-(-P + (P P k2)+{P P ki) 


Plugged into Eq. (El this is matched onto 


(63) 


rjiO^P _ 

■^Aq — 


du6{P+ - u;)—bn^Qfu)r^Y'^%, 


with the four-quark operator 

= 9=5;'23(a;)6+(i..i)7(a;2)(/('‘)'5(Pi(r/)'6+(i..3). 


(64a) 


(64b) 


Here, 5+(a;) acts on everything on its right except /. We include the g"^ in the operators, 
and do not think of it as dvra^, but treat it like the g"^ inside {}Df, as suggested in Ref. ji^ . 
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B. The Result for General Currents 


From the general basis of light-cone operators in Eq. O we define the combinations 
'^^8±(^) ~ g (cn), 

The operators 7^4,io±(<^) are only needed for B —> Xuii^e, bnt not for B —>• Xs'y. Their major 
contribntions arise from the “glnonic” parts 04 ^ 10 ( 0 ;). Employing Eqs. the zero- and 
one-glnon matrix elements of the operators are straightforward to calcnlate, and comparing 
with Eqs. (ESI) and dnn) we can read off their OPE coefficients. 

To write down the light-cone OPE of T“^, we write it as 

= bA^T°-f>^T^b + T"/. ( 66 a) 

The fonr-qnark contribution is given in Eqs. ()64|1 . The expansion of to reads 


( 66 b) 

+ 0{e^). 

For comparison, we wrote the terms in the same order as the corresponding ones in Eq. dSSD- 
Eqs. fj 66 |l represent the light-cone OPE of to (P(£^) and are the starting point for the 
further analysis. 


= ; / da;(5(P+ -o;)<' n^(Po(cu) + 


1 

'K 


n^{0 - V),±{u;) + —n>^{0 + V)s±{u;) 


p2 '^10± 


C. Semileptonic and Radiative Currents 

To continue we consider the currents in Eq. m- For f = u, corresponding to 
B —>■ Xuii^i, we have F“ = 7 "Pl, and the Dirac structure in Eq. (I66a|l becomes F = 
po^Mp/3 = Therefore, with the help of Eq. (IHHll . 

n^F“7'"F^ = (n"n^h,, - [t] + E)“^ n^)YPL, 

h/,F"7''F^ = {n'^n^nu - (7 - ie)'fn^)YPL, (67) 

Similarly, for B —>• Eq. dTUll gives F" = y^^P^. Using = 0 and Eq. dnUl) yields 

= -2n'^(r] + ie)fn,YPR- (68) 

In this case only the terms proportional to in Eq. (I66bfl contribute. 

The appearing Dirac structures are F = 7 “Pl and F = j'^Pr. The discussion in Sec. IIIIBI 
shows that the 7“75 part of and the 7 " part of vanish in the B expectation 
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value, and can thus be dropped. The same is true for the parity even and odd parts of the 
7^4 io_i_(a;), which only come with T = 7 "Pl. From Eqs. we dehne 


(69) 


n^SJcj) = “(u,) + 

The sign of the parity-odd parts in 'R^iQ^ioj) has changed due the minus sign from Pl. 

We parametrize T"^ in analogy to the hadronic tensor in Eq. as 

T°‘^ = J duj6{P+—uj)^—^{ri + ie)^ti — ^{7] — ie)'^t 2 —n‘^°‘n^h 3 + n‘^n^t 4 + n°'n^t^, (70a) 

where the structure functions ti are scalar functions of u and P_. They can be read off from 
Eq. () 66 bfl using Eqs. (inzD and (El. For B —>■ Xuiui we obtain 


= -Ur, 


= -n^ 


fi = - 

^ 2P. 


Ot(^) - ^(o + T)U^) 




fU _ 

Ln - 


P. = 


2 P2 

1 

Ip 


naiO-V)^Au), 


(70b) 


The Pl gives an additional factor 1/2 for the Of (a;) and P2,io(^)^ —1/2 for V^giuj) 

compared to Eq. (lOObjl . For B —> only is nonzero, while the t |_5 vanish identically. 


tl = 2n„ (cn) - — (O - P)?l(^) j , tU = 0- (70c) 

Note that the only operator structures for / = s are Oo, 5 (<^)- In particular, as in Ref. ji^ . 
03 ^ 4 ( 0 ;) do not appear in the QCD light-cone OPE. In both cases they arise only if the QCD 
light-cone operators are expanded into HQET ones. 

Finally, we consider the four-quark contribution T^g . It contains the same Dirac structure 
F = F'^y'^F^. Therefore, from Eq. ()64bfl we dehne 




(71) 


We included a minus sign in the second dehnition, because parity only allows the axial part 
of the left-handed light-quark bilinear to contribute. The four-quark contributions to the 
B —>• Xuii^i structure functions are 


fU _ 

''45,1 — 


1 


2P_ 


n«(Qf + Q^)“(o;), 


fU _ 

''45,4 — 


4P_ 


h„(Qf + Q^)“(o;), 


(72a) 


and for B —> X^y 




(72b) 
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Eqs. d7n|l and (f7^ provide the light-cone OPE of to C>(£^) for B Xuii^e and B —>• Xg'y. 

As a nontrivial check we tested our results against the known local expressions by plug¬ 
ging the shape-function parametrizations of the operators into Eqs. ((70 1 and employing 
their moment expansions. This reproduces the full local result to 0{X) We also 

checked all local 0{X‘) contributions that should be fully contained in our results with the 
expressions in Ref. js^. We hnd agreement for the form factors in the notation of 
Ref. js^. Concerning T^\ we actually disagree with Ref. We believe the contribution 
proportional to 1 /Aq, corresponding to a subleading twist term, should read 


2 (pi -|- 3^2) 
3mbA^ 


(2mb q ■ v){mb 


q -v) + 


2p2mb 

A^ 


We explicitly verihed this by directly computing this term in the local OPE. Ref. js^ misses 
the 2mb in the hrst term. This might have been overlooked so far because is not needed 
in the decay rates for massless leptons, only for r leptons. 


V. DIFFERENTIAL DECAY DISTRIBUTIONS 

We will now use our results from the previous section to derive expressions for various 
differential decay distributions. The decay rates in this section can be used over the entire 
phase space to study arbitrary cuts on kinematic variables. They are valid to 0{X) for small 
hadronic masses sh ~ and to 0{X) for large hadronic masses sh ~ 

For practical purposes this holds provided all shape functions are modelled with correct 
moments up to 0{X). In the resonance region, the expansion necessarily breaks down, and 
one has to integrate the rates over a sufficiently large region to trust the results. 


A. Hadronic Tensor 

It is straightforward to take the B expectation value of Eqs. dZOl) and dH and use the 
parametrizations in the appendix to express the hadronic tensor in terms of the full set of 
shape functions appearing to 0(6“^) at the operator level. However, this includes many higher 
order corrections beyond subleading twist and second local order, as we saw in Sec. IIII Bl 
For phenomenological purposes it is more desirable to reduce the number of shape functions 
as much as possible. 

Since our expansion is accurate to subleading twist and second local order, we can neglect 
all shape functions of twist (!2(A^). In addition, it suffices to keep only those subsub leading 
shape functions which have moments at 0{h?). The hnal parametrizations of the operators 
appearing in the OPE including these simplihcations are given in Eqs. (IA.13|1 . 

Writing the structure functions Wi in Eq. (uni) as 

W(P+, P_) = J duj6iP+ - u;)w,iuj, P_), (73a) 

and using Eqs. (IA.13(1 to take the B expectation value of the in Eqs. (I70bfl . the hadronic 
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tensor for B —>■ Xuii^e becomes 

'^1 = 2 ^-^0 ~ Kq — — (G 5 — (a;), = - ^Fq + Kq — — (G 5 + (a;), 

W2 = ~ “ '^)-^o + -R4] (<^), “ '^)-^o ~ -R4] (i^), 

< = ^ (^4 - 2 (a; - \)Ko - [Rw - AF 4 + 2 A(o^ - A)Fo]^ ( 0 ^). (73b) 

The w“ are regarded as functions of a; and P_. For B —> Xg'j we set P_ = Mb, and Eq. (I7()c|l 
yields 

( 0 ;) = 2 (^Po + Po - ^(^5 - P 5 )) ( 0 ^), w^i^) = 0. (73c) 

Expanding all shape functions into HQET shape functions to subleading twist, the 
vanish and the ^“35 and wl reproduce the expressions for the hadronic tensor in Refs. [23 
li^ . The additional 0{X) contributions have not been computed before. 

Considering the four-quark operators, we define [see Eqs. (IA. 6 fl ] 

{S{“H)b = G{Mk” + --., (S^”M)B = i/|H(n-!))” + ..■ (74) 

The four-quark shape functions are different for charged and neutral P-mesons and f = u,s. 
To avoid having to distinguish between these cases we will not consider them explicitly 
in the following. If desired, they are incorporated by replacing 6 * 5 ( 0 ;) —>• ( 6*5 -|- G'{)(a;), 
H^{ijo) —> (P 5 -|- hI){uj) in all expressions. 


B. Shape-Function Models 

To illustrate our results in the next subsection we employ three models for the shape 
functions, based on the two model functions 

TT r TT fUJ\ 

Pmodi(o;) = ca^o;exp -- j 

n 

Fr.o^2{uj) = c-^o;“'-'e-“-0(o;) 

r(ao) 

The first function is an extension of the one given in Ref. [ 2 ^ and is used in the first model. 
The second one is taken from Ref. and is used for the second and third model. The 
moment expansions of the shape functions are given in Eqs. ()A.14fl . (lA.lhfl . and ()A.16fl . 
Note that the moments are taken with respect to a; — A. 

The leading shape function Fq{uj) is modelled from Pmodi, 2 (i^) by adjusting the parameters 
a, b, c to produce the correct zeroth, first, and second moment. For this purpose we set 
Ti = T 2 = 0 and use Ai = —0.27 GeV^, A 2 = 0.12 GeV^, and = 4.65 GeV, corresponding 
to A = MB—TTih = 0.63 GeV, as our default values. They are inspired by the values obtained 
in Refs. jl| and 13 The third moment of Fo{uj) predicts pi = 0.055 GeV^ in model 1 and 
Pi = 0.084 GeV^ in model 2. The left plot in Fig. HJ shows Fq{u}) in model 1 and 2 . Fo{u}) is 
the same in model 2 and 3. 
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6(uj) + c(l - a) exp 


4 /a;\ 2 


TT 


d{uj), 


(75) 
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FIG. 4: (color online) Shape functions. Model 1, 2, and 3 are dark (violet), medium (orange), and 
light (green). The left plot shows Fb(a;). The right one shows G^{uj) (solid), H^{uj) (dashed), and 
Ri{oj) (dotted). 


The zeroth moments of the subleading shape functions ^ 5 ( 0 ;), and i? 4 (a;) vanish 

to all orders in 1 /mb, because the functions arise from operators containing the derivative 
Therefore, it seems natural to model them by the derivatives F^odi 2 i^)- 
first model we set a = 0 , to ensure that the functions vanish at a; = 0 , and adjust b and c to 
reproduce the correct hrst and second moments. For the second moment of i? 4 (a;) we use 
P 2 = —O.OSGeV^. The second moments of ^ 5 ( 0 ;) and H^{ijj) vanish at C)(Aqqq). We set 
them to (2/3)(0.5 GeV)'^/mb and (0.5 GeV)"‘/mb, respectively. 

In the second model we adjust a, b, c such that the functions have the same hrst, second, 
and third moment as in model 1. Their shape is actually quite sensitive to the value of the 
third moment, which is of CI(Aqqj^). In our third model we adjust the third moments of 
G 5 (a;), and i? 4 (a;) to one half their values in the second model. Therefore, model 2 

and 3 differ only in the subleading shape functions’ third and higher moments, which are 
CI(Aqcd) and higher. The three functions are shown for each model in Fig. |31 Throughout 
the paper we plot model 1, 2, and 3 in dark (violet), medium (orange), and light (green), 
respectively. Notice that ^ 5 ( 0 ;) and H^{ijj) behave roughly oppositely, which means the 
combination (G 5 — H^){uj) is rather large, while (G 5 + is small, as one would expect 

from their hrst moments. 

For modelling purposes we set [see Eqs. (ITT|l and (IA.12f) ] Kq{uj) = [{ijj — \)FQ — R^{oj)/mb, 
where in this case the hrst moment of Fo(a;) is set to zero, and i?io(a;) = — (a; — \)H^{(jj). 

Note that in our approach we regard the local parameters, e.g., Ai and mb, as known 
input parameters. To a hrst approximation, the error due to the uncertainty in their values 
should be treated separately from the error due to the unknown form of the shape functions, 
i.e., their unknown higher order moments. For example, the total rate is very sensitive to mb, 
but basically shape-function independent. Therefore, we do not vary Ai and mb over a large 
range to produce diherent shapes for the functions, because this exaggerates the uncertainty 
for any quantity which is dominated by the local expansion or the hrst few moments of the 
shape functions. Instead, to assess the true sensitivity of a given quantity to the specihc 
form of the shape functions, we look at its variation between diherent shape-function models. 


® For R 4 {uj) this can be seen, for instance, from Eq. (IA.9II . 
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while keeping the zeroth, hrst, and second moments of all shape functions hxed. To do so, 
we use the three models described above. In a more extensive treatment one should include 
more model functions and also scan over generic values for the higher order moments. Of 
course, it is not possible to completely disentangle the two uncertainties, for example, the 
shape-function dependence itself might be different to some extent for different values of Ai 
or rufe. 

Our main interest is in if and by how much a given quantity is influenced by shape 
functions effects and how sensitive it is to the specihc form or higher moments of the shape 
functions. We emphasize that the variations we will see in the plots do not represent total 
uncertainties. They give a measure of the shape-function dependence alone, and are only 
one, in some cases small, part in the total uncertainty. To estimate the latter, one has to 
vary the local parameters as well. To illustrate this, we will vary in the range ±50 MeV 
in a few cases. 


C. Decay Spectra 

We are now ready to assemble the expressions for various decay spectra. Since there 
are no phase space restrictions on our results, we can easily switch to any desired set of 
kinematic variables. We use the notation 

A = Tfg — 2T/£, = Mq — (u, = Zk — (jj = — 2T/£. 

The spectra in all plots are normalized to the partonic rate T" = T^m^ or T^ = T^m^, 
respectively. 


1. Photon Energy Spectrum 

We start by writing down the B —>■ Xgj photon energy spectrum. From Eqs. © and 
(j73cp we have 


dr® / 1 \ 

— = 16rS£“9(i5,) i^F„ + K„ {Mb - 2E,). (76) 

We do not expand the overall E^, because there is no need to do so. Interestingly, Eq. dzni) 
does not contain any subsubleading shape functions. Instead, Fq{uj) and Kq{u) already 
contain all local 0{h?) pieces that are of subsubleading twist. The photon energy spectrum, 
normalized to the partonic rate T® = T^ml, is shown in Fig.|Sl Since at tree level, its support 
lies entirely in the shape-function region, our expansion yields only a small correction of 
0{ls?) to the subleading twist result. 
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FIG. 5: (color online) spectrum. The solid lines show the full result, the dashed lines include 
only the contributions proportional to Fq{u)), and the dotted ones show the subleading twist result. 


2. Spectra in P+ and P- 


With the = Wi{uj, PJ) given in Eqs. ()73bfl . the triple differential decay rate, Eq. (ED, 
becomes 


d3r^ 




= 48K / dw 


S{P+-u;) 

P.-UJ 


e{AAe{p. - A)e{MB - P-) 


X MAMb - P-) (A - P-) + A>2 “(o;, P-) - 2AAA - P_) 


(77) 


X 


M^Mb - P-)W3 (^, P-) + ^>4 (^, P-) + (Mb - P-YwlA, P_) 


The double differential rate, Eq. (HHD, reads 
d2r« r 


dP+dP_ 


= 8r“ / duj6{p+ - uj)e{p_ - uj)e{MB - P-)(P- - ujy 


X 


(78) 


MAMb - P-){< + y^2+ w^3)(^, P-) + P_) + {Mb - P-)X(^, P-) 


To use q± one has to replace P± = Mb — q±- 

The spectrum in the variable P+ is interesting, since it can be directly compared to the 
photon energy spectrum in P ^ to determine the ratio |14fe|/|h)s| Integrating 

Eq. d7H|l over P_, we obtain 

dr“ f r / 1 

^ = KJ du6{P^-u)e{MB-u)lM^^lFo{u) + -Ko{u;) 

2M„ 


Mi ^{—+ 3uj{Mb + oj)) T QMbuj \u.{oj/Mb) ) (Gs — — P 4 — L^){uj) 


+ [-MAM^ - 2u;) + Hu/M b)) {G, + 

— 2Mi^ ( M^{Mb + 5u) + 2u{2Mb + u) ln(a’/iWij)'j [Gg — Hg + Pio — A(P4 + LHiH 


— 3 P 12uMb) + 6Mbu{Mb + u) ln(a;/ilTB) j (Gg + Hg){u) (79) 
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FIG. 6; (color online) P+ spectrum (left) and its partial rate (right). The meaning of the curves 
is the same as in Fig. The vertical line denotes = rnj^/MB, the kinematic limit for b ^ c 
transitions. The horizontal line on the right denotes the value of the total rate to 0{h?‘). 

with [see Eqs. (TOTill . (TOTTl ] 

= —‘2(uj — X)Kq(uj), Gs{u) = —2(a; — A)^Eo(^)) Hq(u) = (o; — A)i?4(a;). (80) 

Notice that the leading term comes indeed with a power M^, as suggested by the subleading 
result, conhrming the leading order result obtained in Ref. js^. The integration over P+ 
to obtain the total rate amounts to dropping the S{P+ — uj) under the inte gra l in Eq. 
Expanding Eq. (|7^ to subleading twist reproduces the results in Refs. |27l. l28|^. 

The P+ spectrum (I70|l normalized to T^ is depicted on the left in Fig. El As for the 
photon energy spectrum, at tree level it has only support in the shape-function region. 
The corrections of the full result (solid) to the subleading twist result (dotted) are thus 
small, although larger then in case of in Fig. El The right plot shows the corresponding 
partial rate, i.e., the spectrum integrated up to P+ < P™* (still normalized to r“). Beyond 
pcut _ pg GeV the curves stay constant. One can see that independently of the used model 
our result indeed approaches the value for the total rate including 0{h?) corrections. 


3. Spectra Containing the Hadronic Invariant Mass 

The hadronic invariant mass sh is useful for the extraction of \Vub\ Using 

SH = P+P-, = {Mb - P+){Mb - P-), dsrf = Mb{P- - P+)dP+dP_, 

and dehning = sh/oj the triple and double differential rates in terms of sh and are 

= 48r“ f - M^{Mb- _ A)e{MB - s^) 


dEidsndq"^ J a;(s^ —a;) 

X Im^{Mb - Suj)\{A-s^fw'^;{uj,s^) + AIw2{uj,Su;) - 2A^{A - s^) (81) 


X 


M^{Mb - S^)W3 (a;, s^) + M‘^w^{uj, s^) -k {Mb - s^fw^{uj, s^) 


^ In comparing our result with Ref. we set the additional P_ cut employed there to zero. 
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FIG. 7: (color online) sh spectrum (left) and partial rate (right). The solid lines show the full 
result, the dashed lines only include the contributions from Fo(a;), and the dotted lines correspond 
to the prescription of Ref. |^. The vertical line denotes sh = The horizontal line on the 

right denotes the value of the total rate to 0{ls?) in units of T". 


and 


ds^dg" 


= 8r“ / duj5[q^ - M^{Mb - s^)]9{s^ - u)e{MB - s^) 


- uf 


00 


X 


M^(Mb - s0l(m“ + «>; + w^)(ui, s„) + M^w^(u}, s„) + (Mg - s„) s„) 


(82) 


Eqs. m and (|H^ can easily be integrated to give d^r“/ dE^dsH and dr “/dsn by dropping 
the 6 fnnction. In this case, the phase space limits yield the limits on the oo integration 


n ^ . , / sh/^ 

0 < -VT- < a; < < ^ , 

Mb \Ash/sh 


for ^/sh < a < Mb, 

for sh/Mb < a < ^/sH < Mb, 


0 < 


Sh 


(83) 




B 


<00 < 


for d‘^T^/dEedsn and dF^/ds/f, respectively. For the latter, upon integration over oo, the 
limits on s/f are 0 < s/f < M^. 

The hadronic invariant mass spectrum dr^/ds^ obtained from Eq. (I82|l is shown on the 
left of Fig. [71 The right plot contains the partial rate for an upper cut sh < The 
solid lines correspond to the full result. The dashed ones show the result keeping only the 
contributions from Eq{oj), in which the complete rrib dependence of the partonic spectrum is 
convoluted with ho (a;). For comparison, the dotted curves give the result from using Eq{oo) 
in the prescription of Ref. ji^®, where an overall ml is excluded from the convolution. The 
full result lies in between the two and neither gives a better approximation than the other. 
One should also note that the solid medium and light (orange and green) curves only differ 
in the third and higher moments, which are (^(Aqcd) and higher, of the subleading shape 
functions. In particular, they share the same medium (orange) dashed and dotted curves. 

Expanding the sh spectrum obtained from Eq. (IH^ to subleading twist reproduces the 
result of Ref. In Fig. IHl we compare our result (solid) for the spectrum and the corre¬ 
sponding partial rate with the subleading twist result (dashed). In addition, the dotted lines 


For consistency this includes only the tree-level results of Ref. |4 
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FIG. 8; (color online) Comparison of sh spectrum (left) and partial rate (right) with subleading 
twist results. The solid lines show our result, the dashed lines the subleading twist result of 
Ref. 241, and the dotted lines the result of Ref. [i^. The vertical line denotes sh = 'mjj. The 
horizontal line on the right denotes the value of the total rate to 0{h?) in units of T! 
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FIG. 9: (color online) Event fractions for sh spectrum. The left plot corresponds to the right one 
of Fig. [3 and the right plot to the right one of Fig. |H1 


show the result of Ref. [2^, which keeps certain factors of and hence contains some 


higher order terms compared to the result of Ref. j2^. The corrections to the subleading 


twist result from our result are more significant then in the or spectrum. For the 
partial rate, in the third model (light, green) they are bigger then the difference between 
the individual models. Our result also has less sensitivity to the form of the shape functions 
then the result to subleading twist. 

The solid lines in the right plot in Fig. |H1 approach the horizontal line, which is again 
a manifestation of the fact that our result contains the total rate to For practical 

purposes, the partial rates are usually translated into event fractions by normalizing them 
to the respective predicted total rate. This introduces an additional error, if the total rate 
is not reproduced correctly. Fig. El shows the event fractions corresponding to the partial 
rates on the right of Figs. [7|and|Hl 


Eq. (IH^ also allows us to obtain the sh spectrum with an additional lower cut on q 


proposed in Ref. 3|. Fig.ITTIshows the spectrum and its partial rate for the cut > 8 GeV^ 
employed by BaBar and Belle Hi- As expected, the cut on high signihcantly reduces 
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FIG. 10: (color online) Phase space cuts. Solid, dashed, and dotted lines are constant sh, Ejj, and 
q'^. The dark (violet) ones are sh = rnjj, Eh = rno, and = {Mb — rnDY'- The medium (orange) 
ones are sh = (1.7GeV)^, Eh = 2.15 GeV, and q^ = 8GeV^, and the light (green) dashed line is 
Eh = 2.45 GeV. 
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FIG. 11: (color online) sh spectrum (left) and partial rate (right) with a cut > 8GeV^. The 
meaning of the curves is the same as in Fig. [71 The vertical lines are sh = w-l) and sh = (1-7 GeV)^. 


the shape-function dependence compared to Fig. | 7 | This comes at the price of having a 
much smaller number of contained b ^ u events. In Ref. 3| the correction from smearing 
the local result with the leading shape function is translated to 100% into an uncertainty 
on the partial rate. Our results can be used to improve on that, and additionally allow one 
to include 0{K) corrections. 

Another possibility is to replace the cut by a cut on the hadronic energy Eh = 
(P+ -|- P-)/2. The rates in terms of Eh are obtained by replacing 5[q^ — M^{Mb — s^^)] —>■ 
5[Eh — (s^^ -f-a;)/2] in Eqs. (IHTll and ([H^ . A cut on Eh < Mb — \/ (g 2 )min produces the same 


upper limit on sh as q^ > (g^)™™. 

The various phase space cuts are depicted in Fig. (TUI The solid dark (violet) and medium 
(orange) lines are sh = sh = (1-7 GeV)^. For (g^)™“ > 8 GeV^ (medium or orange 

dotted) this corresponds to the cut Eh < 2.45 GeV (light or green dashed). The spectrum 
and partial rate for this cut are given by the dashed lines in Fig. [T21 For sh < (1-7 GeV)^ 
the Eh cut signihcantly increases the fraction of 6 —*• u events from about 30% to about 
50% (at tree level), while achieving almost the same shape-function independence. On the 
other hand this cut has less separation power then the q^ cut to reject contamination from 
b ^ c transitions. Lowering it to the intersection of = 8 GeV^ and sh = (1-7 GeV)^ yields 
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FIG. 12: (color online) sh spectrum (left) and partial rate (right) with cuts Eh < 2.15 GeV (solid), 
Eh < 2.45 GeV (dashed), and for comparison > 8GeV^ (dotted). 



FIG. 13; (color online) ruf, dependence of partial rate for sh < ■ Here, dark (violet), medium 

(orange), and light (green) lines correspond to mj, = 4.6, 4.65, 4.7GeV, respectively, while the 
three different models are dotted, dashed, and solid. The left plot contains no additional cut, 
corresponding to Fig. [Tj The right plot corresponds to Fig. where the upper set of curves is for 
Eh < 2.45 GeV, the middle set is for Eh < 2.15 GeV, and the lower set is for > 8GeV^. 


Eh < 2.15 GeV (medium or orange dashed). This still retains a larger fraction of the signal, 
and at the same time cuts out a somewhat larger portion of the 6 —> c phase space, which 
should in principle provide an equal or better suppression of contamination from b ^ c. 

The Sh spectrum and partial rate for this cut are shown by the solid lines of Fig. IT^ 
For comparison, the dotted lines show the result for the cut, i.e., they are identical to 
the solid lines of Fig. The Eh cut has basically the same shape-function independence, 
but retains additional 5% of signal events and should provide an equal or better 6 —>• c 
separation. Of course, eventually this depends on the experimental resolution. We conclude 
that a combined analysis of sh and Eh provides a viable alternative for measuring iKiftl, 
with a potentially higher accuracy than the analogous measurement of sh and 

Fig. IT^ shows the variation of the partial rate for sh < for the various cuts when 
changing rrib by ±50 MeV, where curves of the same color correspond to the same value of 
mb- It shows, that the uncertainty in mb is clearly a separate, in these cases much bigger, 
effect then the sensitivity to the specihc form of the shape functions. The two effects should 
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therefore be distinguished and treated as separate uncertainties, as argued at the end of 
Sec.IVU 


Hadronic Energy Spectrum 


A hxed hadronic energy corresponds to a line with slope —1 in the P± plane, see Fig. ITIIl 
The hadronic energy spectrum thus receives contributions from both local and twist phase 
space regions, and is therefore interesting to study in its own right. In addition, considering 
the single differential spectra, a cut on < rriD (dark or violet dashed in Fig. ITn|l retains 


substantially more signal events than the equivalent cut on 


dotted). Both methods have been suggested to determine |I4fe| |45Ll4(l l4 

Changing variables from to Eh in Eq. (jS^ and integrating over sh, we obtain 



(dark or violet 


dF“ 


= 64F“ y dueiEn - uj)e{MB - 2 Eh + uj){Eh - ojf 

M^{Mb - 2Eh + a;)« + < + 2Eh - u;) + 2Eh - oo) (§4) 


X 


+ (Af^ — 2Eh + 2Eh — Co’) 


Using the fact that 0 < cu, the phase space limits yield the integration limits 

0 < O’ < Eh for 0 < Eh < Mb/2, 

2Eh — M ^ to ^ Eh for Mb/2 ^ Eh "E Mb- 

The Eh spectrum and the partial rate obtained from integrating it up to Eh < are 
shown in Fig. Our result (solid) matches the local result (thin black) over a wide range 
of energies and smooths it out near the partonic phase space boundaries. The partonic 
boundary however does not he in the shape-function region, and hence, the shape-function 
dependence in the dropoff at Eh = mb/2 -|- A, i.e., the differences between the three models, 
are very mild, for example compared to the sh spectrum or the lepton energy spectrum (see 
below). In the Eh spectrum the prescription of Ref. seems to give a better approximation 
than convoluting the full mb dependence of the partonic spectrum. A cut on Eh < thb alone 
keeps 21% of the b —>■ u signal (normalized to the partonic rate), which is 50% more than 
the cut on > (m^ — m^i)^, which keeps 14%. This relative increase should not be changed 
much by radiative corrections. At the same time the spectrum and partial rate in this 
region are completely shape-function independent. Hence, measuring the hadronic energy 
spectrum alone to extract |I4b| seems worth pursuing, too. 

The spectrum and partial rate with an additional cut sh < m|, (solid) and sh < 
(1.7 GeV)^ (dashed) are given in Fig. HSl A slight increase in the sh cut allows one to 
substantially raise the Eh cut while still keeping the partial rate practically shape-function 
independent. Ideally, if the cut sh < m|, would remove all charm background, the Eh cut 
could be raised up to Eh < 2.7 GeV, which would yield a partial rate around 70%. The mb 
dependence of the partial rate with these cuts is shown in Fig. The uncertainty in mb is 
again important, as one would expect, but distinct from the shape-function uncertainty. 
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FIG. 14: (color online) Eh spectrum (left) and partial rate (right). The solid lines show the full 
result, the dashed lines only include the contributions from Fq{u;), and the dotted lines correspond 
to the prescription of Ref. |^. The thin black curve is the local 0{h?) result. The vertical line 
denotes Eh = mo ■ 



FIG. 15: (color online) Hadronic energy spectrum (left) and partial rate (right) with a cut sh < w-f) 
(solid), Sh < (l.TGeV)^ (dashed), and no cut (dotted). The thin black curve is the local 0{h?‘) 
result, and the vertical lines denote Eh = m^). Eh = 2.15 GeV, and Eh = 2.45 GeV. 



FIG. 16: (color online) mb dependence of partial rate for Eh < E^^ with a cut sh < 'm'jj (left) 
and Sh < (1-7GeV)^ (right). The meaning of the curves is the same as in Fier. II,‘IL 
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FIG. 17: (color online) spectrum without cut on Eg (left) and with Ei > 2GeV (right). In 
each case, the smooth black curve is the partonic result and the one with the edge is the local 
result. 




FIG. 18: (color online) spectrum (left) and partial rate (right) with Ei > 2 GeV. The solid 
lines show the full result, the dashed lines only include the contributions from Eq{uj), and the 
dotted lines correspond to the prescription of Ref. Q . The thin black curve is the local result to 
O(A^), and the vertical line denotes = mj^. 


5. q^-Ei Spectrum 

The q^-Ei spectrum is also of great interest for measuring |14fe|- Integrating Eq. (IHT| over 
sh, and defining = q'^/M^, we obtain 


d2r« 

dEfdg" 


= 48r) 


do; 


M^{M^ - q^ 


-9{A^)9{2Ee - qu)9{quj)lq‘^ {q^j - 2Eifwl{ijO,MB - quj 


+ A^W2 (w, Mb - quj) - 2A^{q^ - 2Ee) q^w^{u, Mb - q^j) 
+ M^wKlj, Mb - quj) + {qujfw^{uj, Mb - q^ 


( 86 ) 
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The resulting integration limits on oj are 


0 < w < 


Mb — ‘2E£ 

Mb - qV{‘^Ee) 


for 0 < 2Eg < < Mb, 

for 0 < \/q^ < 2Ei < Mb- 


(87) 


The maximally allowed hadronic invariant mass for given and E^ dehnes the variable 


= {Mb- 


2E, 


{Mb — 2E(). 

/ 


Requiring is an efficient way to discriminate the b ^ c background 


and 


has been implemented by BaBar p. The distribution in provides a quite nontrivial 
application for our expansion, because fixed contains contributions from large which 
should be shape-function independent ji^, as well as large lepton energies, which are shape- 
function sensitive. The spectrum is obtained from Eq. ()8f)|l by changing variables from 
g^ to and integrating over the lepton energy. It is depicted in Fig. El with and without 
an additional cut on > 2 GeV. Above ~ 5 GeV^ the spectrum is determined by 
the local result, which extends to s™*™ = M^, where it goes to zero. Below that our result 
smooths out the local spectrum. The shape-function sensitivity in the lower part of the 
spectrum is somewhat larger than for E^, but still much smaller then for sh or Ei- With 
a cut Ei > 2 GeV the maximum value of is Mb{Mb — 4 GeV) = 6.76 GeV^. Although 
this cut removes a large fraction of the local OPE part of the included phase space, it is still 
low enough that the increase in the shape-function sensitivity is insignihcant. 

For comparison, Fig.fT^shows the spectrum and partial rate with a cut Ei > 2 GeV, where 
the dashed lines only include the contributions from To (a;), and the dotted lines implement 
the prescription of Ref. |4^ . For the partial rate our result yields a sizable correction to the 
latter and also to the local result. The variation between the different models is negligible. 


6. Lepton Energy Spectrum 

Finally, we come to the lepton energy spectrum. Integrating Eq. (ED over P4. and P_ or 
Eq. (|HT)D over g^ we hnd 

dr“ f f 

^ = 4r^0(Ei) I du;e{A^)Mj4E^{M^ + 2A^)(Po - Ko){u;) 

+ 12M^A^ (2p£-|- ln{Ai^/M^)'jK q{uj) 

- 6(2Ei{Ei + 2 A^ To;) - {M^AI ln(A^/M^)}') {G, - H,){u;) 

- 6 M^A^{A^ ln(A^/M^)}'(G5 + H,){u;) 

- 12A^{2Ei - {M^A^ln{AjM^)y){R, + L,){u;) 

- 3Al{M^ ln{AjM^)Y{Gs - Hs){u;) 

- 3 ^ ( 4 P, - ln(A JMJ}") (Gg + Hs){u;) 

+ 6A^{M^A^ln(A^/M^)}"[Pio - A(P4 + L3)]{u;)], 
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FIG. 19: (color online) Lepton energy spectrum spectrum. The solid lines show the result from 
Eq. (IHHl) . and the dashed ones that of Ref. |i^. The upper black line shows the partonic spectrum, 
and the lower one the nonsingular part of the local spectrum to 0{h?‘). 



FIG. 20: (color online) Lepton energy spectrum (left) and partial rate (right). The solid lines 
show the full result, the dashed lines the subleading twist resultj_and the dotted lines the result of 
Ref. 13 • The vertical line denotes the BaBar cut Ei = 2 GeV 


where ^ 3 ( 0 ;), Hs{uj), and L^{ui) are given in Eq. (IHIHl above, and we use the notation 




/M - /(o) 


00 


{/Mr=2 


f{uj) - /(O) -a;/'(0) 


00^ 


(89) 


Expanding Eq. (j 88 |l to subleading twist reproduces the result in Ref. [27[. It also agrees with 
Ref. [ 2 ^. However, for some reason, the authors there divide their result by an additional 
factor Mb — OJ and subtract a compensating term (a; — A)Fo(a;). 

Eq. and the result derived in Ref. j26| agree when both are expanded to subleading 
twist. However, we cannot expect them to be identical, since the light-cone directions in 


Ref. Hi and in the present case are different. The modihed expansion thus retains different 
higher order twist corrections in each case. The direct computation in Ref. 0 yields a 
much more concise result, because it uses a light-cone direction natural to the lepton energy 
spectrum. On the other hand one could say that the present choice includes more higher 
order twist corrections, since it keeps them already at the level of the triple differential 
rate. However, to make this statement precise one would need to compute all subsubleading 
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twist contributions. The two results are shown in Fig. Below Ei ~ 1.8 GeV they are 
dominated by the local spectrum. Above Ee ~ 2.2 GeV they become quite different, which 
indicates that the higher order kinematic twist corrections are important. 

In Fig. 1201 we compare our result (solid) with the subleading twist result (dashed and 
dotted). We note two things. First, the spectrum on the left shows that the twist expansion 
breaks down around ~ 2.1 GeV, and cannot be trusted for smaller energies. Secondly, we 
see again that the higher order twist corrections are important. Interestingly, the variations 
between the different models are significantly larger for the subleading twist result than for 
our result. We already observed a similar but smaller effect in the sh spectrum. This sug¬ 
gests that the endpoint spectrum is to a large extent determined by higher order kinematic 
twist corrections. If this is the case, it would be promising to the \Vub\ extraction from the 
lepton energy endpoint spectrum, where the shape-function dependence is a limiting factor 
in the achievable accuracy. To conhrm this one certainly needs to consider a wider spread 
of model functions. One could also compute the true subsubleading twist corrections to see 
if they have an equally large effect on the spectrum or not. 


VI. CONCLUSIONS 

We have studied the Aqqd/Mb corrections to inclusive 5-meson decays, with attention 
to the radiative decay B —> and the semileptonic decay B —>■ Xuiu^. 

Usually, the twist expansion is valid in the shape-function region and the local expansion 
in the rest of phase space. Following Ref. |2^ we used a modification of the twist expansion 
which avoids the restriction to the shape-function region and yields an expansion applicable 
over the full phase space, except for the resonance region. This effectively provides a smooth 
interpolation between the otherwise separate regimes of local and twist expansion. So far, 
we only worked at tree level. It would certainly be interesting to see how much of our 
approach can be carried over to include radiative corrections, at least for the contributions 
proportional to the leading shape function, since the corrections to the leading twist 
result are known |5^. To extend the matching calculation presented here to order 
one needs to study the renormalization of the relevant light-cone operators in QGD. The 
renormalization properties of the shape functions will be different when they are defined via 
QGD rather then HQET operators. Similarly, the a* corrections to the relations between 
shape-function moments and the standard HQET parameters will differ between QGD and 
HQET shape functions. 

We performed the expansion directly in QGD light-cone operators and gave a discussion 
of the general operator basis appearing at tree level, including its parametrization in terms 
of QGD shape functions. We used reparametrization invariance under rotations of the light- 
cone direction to reduce the number of independent shape functions, and showed that the 
different lepton energy spectra obtained in Refs. |26L 1271 128|| are in agreement. 

The results for the various decay distributions are presented in Sec. 0 The photon 
spectrum for B —> X^y is given in Eq. dZSl). The triple differential decay rate for B Xui-Vi 
is given in terms of 5^, P+, and P_ in Eq. (|77|) and in terms of P^, Sh, and in Eq. (ED. 
Except for the resonance region, the decay rates are valid over the entire phase space, being 
exact to order Aqcb/Mb in the region of hadronic masses sh ^ 0{AqcbMb), and to order 
Aqqd/Mb away from it. 

Employing different shape-function models our results allow to quantify the impact of 
shape-function effects on decay distributions and partial rates for any desired kinematic 
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cuts. We stress again, that the shape-function models we used only differ in the third and 
higher moments of the shape functions, but still give quite different shapes. The observed 
variations in the results provide a direct measure of the true shape-function sensitivity of a 
quantity, i.e., its sensitivity to the unknown specific form or higher order moments of the 
shape functions. For the total uncertainty, one has to additionally vary the local parameters, 
most notably Ai and m^, as well, the effect of which should be regarded separately. 

An application is to study the transition between the local and twist expansion. The 
primary example is the lepton energy spectrum, and we saw that the usual twist expansion 
cannot be trusted below Ei ~ 2.1 GeV. 

We are free to choose any kinematic variables, and discussed several examples of interest. 
In particular, we can study decay rates which for given values of the kinematic variables 
receive contributions from the phase space regions of both local and twist expansion, such 
as the hadronic energy, or the variable used in the q^-E^ analysis. The hadronic energy 
spectrum has not received much attention so far. We point out that, with or without an 
additional cut on sh, it represents a viable alternative to the existing SH-q^ analyses to 
extract | 14 b|- 
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APPENDIX: PARAMETRIZATION OF LIGHT-CONE OPERATORS 


Here we collect the results related to the parametrization the light-cone operators from 
the body of the paper and add some further details. First, the number of independent 
Lorentz structures in the operator basis in Eq. O can be reduced employing the relations 




= 0 , 

= -a;C>4^(a;), (A.l) 
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They 


(A.2) 


which hold for any Dirac structure T. Note that there is no relation for 
also imply 

(uj) = = -[a;^C>o(n;)]', 

n^n^Ol^''{uj) = n^nyOl^''{uj) = a;^C>o(a;), 
n^nvO^^''{uj) = n^nuO[l^''{uj) = 0 . 

For completeness we repeat the parametrization of the leading operator, Eq. (HOD, 

+ iFoM(n - n)“, (A.3) 

and the 0{e) operators, Eqs. PH|) . 

= -{F, - AFo)'(o;)n“ - {K, - XKoy{u;){n - n)“, 

(h±a/.C>r(^))i? = 

{v,Oy>^{u;))B = (Fs - AFo)(a;)n“ + {K, - XKo){u;){n - n)“, (A.4) 

{v±at,0^''{uj))B = Lsiu), 

The set of all nonzero matrix elements for the 0{e‘^) operators is given in Eq. (j47j) . We need 

{u))b = G5,8(a;)n" + M 5 ^si^){n - n)", 

= F5,8(a;)(n - n)“ + A^5,8(a;)n“, 

i^^^B = {Lg - XLs)^^), 

= (Fio - AF 4 )(a;). 

The four-quark operators, dehned in Eq. (HU give rise to the shape functions 

(Q(“M)n = Gi{u)v^ + M( H(n - t;)“, 

(Q^“(o;))b = Hy{uj){n - v)° + Ni{uj)v°. 


(A.5) 


They depend on the flavor of the hnal-state quark, and are thus different for B 
B —>■ Xuii^e- In addition, they also differ for charged and neutral F-mesons. 
The RPI constraints in Eqs. m require 


(A.6) 

Xg'j and 


(A.7) 


(A.8) 


L[{uj) = 2Ko{uj), — 2Fi(a;), 

^ 2 ( 0 ;) = -2[(a; - A)Fo + Fi](a;), M'(a;) = -2[(o^ - X)Kg + K,]'{u:) - 2iFoM, 

where the functions (^2(1^), L2{uj), and M2 (a;) are dehned by 

{v±f^i'G>2^^''\^))B = -IG2 {uj)v°‘ - |M2(a;)(n - n)“, 
{V±a(^V,)Of>^‘'\uj))B = \{L2 - XL^)"{u). 

They also require 

{VL^.uOy^\uj))B = -2[{uj - X)Fo + Fs]{u;)v^ - [2{u - X)Kg + 27^3 + Ly\{uj){n - n)“, 
{v,v^auOt^''{^))B = -(F 3 - AFo)(o;), 
yie^^yVy^''{uj))B = Ri{uj){n - n)" 0n“, 

{v,ie^o..Vy^\u2))B = 0 . 


(A.9) 
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In the remainder of this appendix we use A" to denote Aqqq divided by an appropriate 
power of nib. The functions and Lg{uj) are suppressed by 1/mb. Hence, 

they are effectively twist Cl(A^) and can be neglected. The HQET equations of motion imply 

F,{u) = OiA^), 

Gsiu) = -2{u; - X)^Fo{uj) + 0{A^), (A.IO) 

= (uj — A)i?4(a’) + Cl(A^), 

where the neglected terms are of higher twist order, and the relations for G^{uj) and 
follow from the form of their n**^ moments. In addition, we may neglect also all twist 0{A^) 
shape functions without moments of local 0{A^). The RPI relations (IA.7jl and (IA.9(1 and 
the HQET equations of motions then imply 

FiM = 0(AQ, ujF[{uj) + FsM = 

K,{u;) = 0{A^), K^{u;) = 0{A^), (A.ll) 

= 2Kq(^uj), L^iuj') = — 2{ijj — A)A'o(ctj) + Cl(A^). 

The neglected terms are now only of higher order in the local power counting. The relation 
for Li(a;) is exact. There is no formal relation hxing i?io(a;), but we may model it as 

=-{oj - X)Hb{uj) + 0{A^) (A.12) 

which correctly reproduces its hrst two moments, see below. 

Putting everything together by employing the above shape-function dehnitions and 
Eqs. fjA.l|l . (IA.10|1 . and (IA.11|1 . the parametrization of the operators in the OPE are 

{Oq{u))b = Eo(a;)n“ KQ{uj){n - u)", 

((^5±M)b = G5(o;K + --- , 

= H,{u;){n - vT + ■ ■ ■ , 

= -2(a; - XfFo{uj)v^ + • • • , (A.13) 

(^8±(^))s = (o; - X)Ri{uj){n - vY ^ -, 

{nTA^))B = \[R,- 2 ( 0 ; - x)Kg]{u;YT + ■ ■ •, 

{KoA^))b = |[^io - XR, + 2A(a; - X)Kg]{u:)tY + ''' • 

At last, we look at the moment expansions of the shape functions. For Fq{uj) and Kq{ui) 
they were given in Eqs. (HBD, 

FYu) = 5{u; - A) - - A) - - X) 

ZTTlfj O 

/ \ 2Ao Po/B^bj-t/ \\ I Po j-n/ \\ I 

7 - 0 (a; - A) + -—b (a; - A) H-. 

omfe omfe 


-^r(o;-A) + ---, 


(A.14) 
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For the shape functions arising from the 0{e) operators we have 


F.iu) = - A) + OiA^)5'iu - A) - ^5"(a; - A) + ■ • • , 

Zuit, io 


FJoj) = -^6(uj - A) + 0(A^)S'(uj - A) + 0(A^)6"(uj - A) + • • • , 

2mb 


K,,,{u;) = -^6 '{u;-X) + ---, 

0772 ^ 


(A.15) 



Finally, the shape functions arising from the 0{e^) operators we need obey the expansion 


— ——(Ai + Ti/mb)6'{uj — A) + 0{A^)5" {lo — A) + ■ ■ • , 

O 

= -(A 2 + T2/mb)5\uj - A) + 0{A^)5'\uj - A) + • • • , 
G%{ui) = ^(-^i + — A) H— ^5'{uj — A) + ■ • • , 


(A.16) 


— (A 2 + T2/Tnb)5(u} — A) + P2A {uj — A) + ■ ■ • , 


— -Z —— A) + • ■ ■ , 

Srub 

Rwiuj) = -(A 2 + T2/mb)6{u - A) + 0{A^)6\uj - A) + ■ • • . 


All moments of and N^^si^) are 0{A^) and higher. 
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